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Surface-electrode (SE) rf traps are a promising approach to manufacturing complex ion-trap net- 
works suitable for large-scale quantum information processing. In this paper we present analytical 
methods for modeling SE traps in the gapless plane approximation, and apply these methods to two 
particular classes of SE traps. For the SE ring trap we derive analytical expressions for the trap ge- 
ometry and strength, and also calculate the depth in the absence of control fields. For translationally 
symmetric multipole configurations (analogs of the linear Paul trap), we derive analytical expres- 
sions for electrode geometry and strength. Further, we provide arbitrarily good approximations of 
the trap depth in the absence of static fields and identify the requirements for obtaining maximal 
depth. Lastly, we show that the depth of SE multipoles can be greatly influenced by control fields. 
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Several approaches to large-scale quantum informa- 
tion processing (QIP) with trapped ions using networks 
of interconnected rf traps have been proposed [1, 2, 3]. 
Surface-electrode (SE) ion traps, where the trap elec- 
trodes are located in a single plane as illustrated in Fig. 1, 
offer a promising approach to constructing such trap net- 
works [4, 5, 6]. SE traps are well suited for microfabri- 
cation, and motional heating rates compatible with QIP 
have been demonstrated [5, 6], an important fact since 
motional heating is currently a limiting factor for ion- 
trap miniaturization [7, 8]. 

Except for analytical results for some high symmetry 
configurations [9] , modeling of SE traps has so far been 
based mostly on numerical techniques. While this ap- 
proach is feasible for modeling with the purpose of char- 
acterizing the trap provided by a given electrode struc- 
ture, it is inconvenient for modeling with the purpose of 
designing traps to best meet design criteria. This is un- 
like the situation in traditional rf traps such as the linear 
Paul trap [10], where the electrodes surround the ion, 
and where symmetry considerations together with scal- 
ing relations often provide sufficient background for trap 
design. 

In this paper, we present two methods for the design 
of SE traps. Firstly, we demonstrate that a relatively 
unknown technique makes it possible to analytically cal- 
culate the field resulting from an arbitrary configuration 
of surface electrodes in the "gapless plane" approxima- 
tion where the electrodes are assumed to gaplessly cover 
an infinite plane. Secondly, we show that the important 
case of multipole guides with translationally symmetric 
electrode configurations (analogs to the linear Paul trap), 
can be given a simple geometrical description through a 
conformal map. Together with recent results on the fun- 
damental constraints on intersections in rf trap networks 
[11], the results presented here can serve as guidelines in 
the design and modeling of SE ion trap networks. 
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FIG. 1: (Color online) Cut near the trap center (oval) of 
a surface-electrode (SE) trap, showing the rf and control 
(DC) electrodes on an insulating substrate. Design, including 
the aspect ratio of the gaps between electrodes, corresponds 
roughly to that of the trap described in Ref. [5], for an ion 
surface distance of 40 /xm. The electrodes extend beyond the 
illustrated area in all directions. 



The paper is structured as follows: Sec. I introduces 
the fundamental task of ion-trap modeling. In Sec. II, we 
argue that SE traps can be modeled as an infinite plane 
gaplessly covered by electrodes and describe an analyti- 
cal method for calculating the field of arbitrarily shaped 
surface electrodes in this approximation. In Sec. Ill we 
apply this method to the analysis of the SE ring trap. 
Sections IV and V are devoted to parametrizing and an- 
alyzing translationally symmetric SE configurations. 



I. RF ION TRAPS FOR QIP 

In this section we introduce the basic physics of rf traps 
and the role of electrostatic modeling in the characteri- 
zation of such traps. 

For the typical dimensions (/im) and rf frequencies 
(MHz) of the rf traps typically used for QIP, the rf field 
can be adequately treated as quasistatic [12]. Follow- 
ing Ref. [13] we introduce the electrode basis function 
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0^(r) for an electrode 7 as the unique function that 
takes the value 1 (0) for r on the surface of electrode 
7 (all other electrodes) and fulfills the Laplace condition 
V^Q^{r) = 0. In terms of the electrode basis functions, 
the spatial potential obtained by applying a voltage to 
electrode 7 while keeping all other electrodes grounded 
is Vj{r) = VyQj{r). Neglecting any nonlinear electric 
properties of the trap materials, the total field ETotif) 
in the trap can be expressed as 

ETot{r,t) = - {Vricos{nt) + Vc,rf) Verf(r) 

-J2v,,iVQ,,i{r), 

i 

where Q. is the rf frequency, Kf is the peak rf voltage, V^Ti 
is the bias applied to the rf electrodes, and Vc^i is the volt- 
age applied to control electrode i. In most applications 
the applied voltages will be modified quasistatically to 
manipulate the trap, but we shall not be concerned with 
this in what follows. To maintain physical dimensions, we 
will express most results in terms of the spatial potential 
Ky(r) and the corresponding field Ej{r) = — VV^(r). 
Note that with this convention E^i is the rf field ampli- 
tude rather than the instantaneous rf field. 

In an rf trap, the rf period 27t/Q is made to be the 
shortest timescale of the system, so that in the adiabatic 
approximation the effect of the rf field on an ion of charge 
Q and mass M is described by the ponderomotive poten- 
tial 



U,{r) ^ 



(1) 



corresponding to the kinetic energy of the micromotion 
induced by the rf field [14, 15]. Including the effects of the 
quasistatic control fields, the effective external potential 
experienced by a trapped ion is 



Ues{r) = Up{r) + Q 



K,rf(r)+^K,i(r) 



(2) 



Although trapping at the local minima of U^s where Up 
is non-zero is possible, it is usually avoided as the micro- 
motion associated with a non-zero rf field amplitude can 
have a range of detrimental effects. In the following, we 
will only consider trapping at points where Erf{r) = 0. 



II. FIELD OF ARBITRARY SURFACE 
ELECTRODES 

In the following, we argue that it is in many cases rea- 
sonable to model an SE trap as an infinite plane gaplessly 
covered by electrodes. Further, we show that in this ap- 
proximation the field of an arbitrary electrode can be 
calculated analytically. Lastly, we consider the field in 
terms of the Fourier components of the electrode shapes. 



A. Gapless plane approximation 

As introduced above, an SE trap has all electrodes lo- 
cated in a single plane, typically on the surface of an 
insulating substrate as illustrated in Fig. 1 [4, 5, 6]. 

We will model the SE configuration in the gapless plane 
approximation, where the gaps between individual elec- 
trodes are assumed to be infinitely small and the elec- 
trodes are assumed to fully c;over an infinite plane with no 
other conductors above the plane. It is clear that with d 
denoting the ion-electrode distance, necessary conditions 
for the gapless plane approximation to be valid include 
that gaps between electrodes are much smaller than d, 
that the extend of the trap is much larger than d, and 
that the distance from ion to other conducting surfaces 
is much larger than d. 

These requirements are often met in SE traps. Firstly, 
although small electrode gaps lead to high field gradi- 
ents, SE traps are usually constructed with small inter- 
electrode gaps since any exposed surface of the substrate 
would be prone to pick up static charges disturbing the 
operation of the trap. The condition on the absence of 
other conductors in the trap region imposes constraints 
on the experimental apparatus, but will in many cases 
be met. The corrections due to nearby conductors will 
be first order in the ratio of the total electrode extend to 
the distance to the conductors. For larger trap arrays, 
this could be a significant effect and corrections should 
ideally be made. 

In the gapless plane approximation, the electrostatic 
modeling for a patch electrode 7, covering a region Aj 
of the electrode plane reduces to identifying the unique 
potential t^(r) so that V^Ky(r) = and 



V-y{x,y,0) = 



for (x, y) e 
otherwise. 



(3) 



where we have chosen the x-y plane to coincide with the 
electrode plane. 



B. Biot-Savart law 

As shown in Ref. [16], the electrical potential Vy{r) 
satisfying the boundary condition of Eq. (3) can be very 
elegantly expressed in terms of the sohd angle QA^{r) 
spanned by A^ as seen from r. 



2tt 



(4) 



Taking the electrode plane to be the x-y plane, and as- 
suming 2; > 0, we have that 



J A 



I A |r-r'| 
The corresponding field takes the form [16] 
„ , , r dr' x(r- r') 

JdA^ |r-r'| 



(5) 



where the integral is counterclockwise (as seen from 
above) along the edge dAj of A^. Note that the inte- 
gral is a Biot-Savart integral, so that is proportional 
to the magnetic field that would be observed if a current 
was run along a wire following the edge of . This points 
to a close connection to work on "atom chips" using mi- 
crofabricated wires [17, 18, 19], but it should be noted 
that this system is fundamentally different since magne- 
tostatic traps can trap at nonzero field minima without 
adverse effects. Numerical integration of Eq. (5) can ben- 
efit from the results of Ref . [20] . 



C. Fourier description of SE fields 

For some applications, e.g. trap arrays [21], it is useful 
to describe the surface electrodes in terms of their Fourier 
transforms [22] . This approach also provides insight into 
the z dependence of SE fields in general. 

By employing the Dirichlet Green's function also used 
in Ref. [16], we find that in terms of fc = k^x + kyij 
the Fourier transform V{k,z) = J e~'^''''^V{r) dx dy of 
V(x,y,z) with respect to the x and y coordinates for a 
constant z value is given in terms of the surface potential 
Vsurf{x, y) = V{x, y, 0) by 



V{k,z) = 



-ihr 



1 



\z-r\ 



2t7 1 7- — 7"' 



3 V;urf(a;', y') dx dy dx dy' 
(6) 




FIG. 2; (Color online) Strongest possible SE ring-trap con- 
figuration. Upper part: Cut plane through the rotational 
symmetry [z) axis, together with the trap center (triangle) at 
a distance d above the electrode plane, cut of ring electrode 
(rectangles below the x axis), and field lines. Lower part: Top 
view of the electrode {x-y) plane, which in the gapless plane 
model is assumed to be grounded everywhere but on the ring 
electrode with inner (outer) radius Ri (R2)- 



where 14urf(fc) = / e~^^'^Vsm:i{x,y) dx dy is the Fourier 
transform of t^urf- Equation (6) implies that the Fourier 
components of 1^ at a given z value fall off exponentially 
with z on the length scale of the electrode dimensions. 

For our purpose, Ksurt will be equal to V^on some region 
A and zero elsewhere. In this case, V^urf can be calculated 
by noting that since e"''" '' is equal to the z component 
of the curl of the vector field i^z x e~*'° '"A;) and the 
vector triple product a x b • c is invariant under cyclic 
permutations, the Fourier transform can be converted to 
the path integral 

V;urf(fc) ^V^z-(k e-^'^-^k X dr, 

JdA 

where the integral is counterclockwise along the edge 
(OA) of A. 



III. RING TRAP 

To illustrate the application of the electrostatic meth- 
ods presented in the previous section to rf trap analy- 
sis, we will, as an example, consider the SE version of a 
quadrupole ring trap [10]. We will show below that this 
can be implemented as a single ring-shaped rf electrode 
with inner and outer radii Ri and R2, as illustrated by 
Fig. 2. 



The field of the ring electrode at potential V can be 
calculated as the superposition of the fields of a disk of 
radius i?2 at potential V and a disk of radius i?i at po- 
tential —V, so that the combined field is -E = Ej^^ — Ej^^ , 
where Eji is the field of a disk of radius R at potential 
V. By symmetry, the electric field on the z axis must be 
aligned along the axis, and by Eq. (4) we find that 



Er{zz) = V 



i?2 



(i?2 + 



n3/2 



(7) 



See, e.g., [23] for the off-axis case. 

For certain values of Ri and R2, we find that E^izz) 
vanishes for some value of z. When rf is applied to the 
ring electrode, such zeros correspond to points of zero rf 
field amplitude, which provide an ideal trapping point as 
discussed in Sec. I. 

Rather than trying to determine the position of any 
field zeros for given dimensions of the ring electrode, 
we will consider the constructive design problem of 
parametrizing all values of i?i and R2 so that a field 
zero is located at a chosen distance d above the elec- 
trode plane. As an intermediate step to achieving this. 
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following we consider the strength and the depth of the 
trap, as illustrated in Fig. 3. 

The strength of an rf quadrupole trap is usually de- 
scribed by the dimensionless parameter defined as [14] 



9z = 9o 



2 



qo = 



(9) 



where go has the same scaling properties as Qz- In the 
adiabatic limit where the ponderomotive approximation 

is valid, the secular oscillation frequency oJz in the z di- 
rection is seen from Eq. (2) to be g^Q/v^. 

For the ring trap the rotational symmetry and the 
Laplace condition imply that the lowest-order multipole 
contributing at the field zero must be a quadrupole term 
of the form E{r) = -af^ {2zz - xx - yy) + 0{r'^). The 

(2) (2) 

quadrupole strength Uz can be calculated as az = 
— I {dEz/dz). Evaluating the derivative of Eji^ — Efi^ 
according to Eq. (7), we can write the resulting strength 

rf^/F explicitly in terms of the parameter 6. 



, (2) 

(Iz = (lo\az 



qz^qo \ [sin(56') - sin(6')] , 



(10) 



FIG. 3: Properties of the SE ring trap as a function of the 
parameter 6 appearing in Eq. (8) for the example parameter 
values given in Eq. (11). Us (top) is the value of UcS at the 
saddle position in the absence of any control fields. uJz = 
Qz^l/VS (middle) is the adiabatic approximation of Wz, valid 
when Qz <C 1. Ri,2 (bottom) are the inner and outer radii 
of the ring electrode. The strongest trap is obtained for 9 « 
0.275, where Cjz ^ 2v x 8.17MHz, Us « 118 meV, and Ri,2 ~ 
{68, 338}/im. Sealing of qz, Us, and i?i,2 are given by Eqs. (9), 
(12), and d, respectively. 



we parametrize Ri as dtan(0j) so that 
Ez{dz) = ^(cos((/.i) [1 - cos2(<^i)] 

- C0S{(I)2) [l - COS^(02)] )• 

For Ez{dz) to vanish, cos((/)i) and cos(02) must be roots 
of the polynomial u(l — u^) — a for some common con- 
stant a. This polynomial equation is in the Vieta stan- 
dard form [24], and the solutions of interest to us can 
be parametrized as cos((/)i,2) = 2sin(7r/6^0) /VS, cor- 
responding to R values of 



Ri,2 = d\j - sm 



< 61 < f. 



(8) 



This, in particular, implies that the smallest outer radius 
of the ring electrode is i?2 = cJ-\/2. 



which reaches a max;imum value of Qz « 0.473 go when 
cos^(0) = (25 + \/145)/40. For typical parameters of 



fl = 2-!TX 100 MHz, 
M= lOAMU, 



= lOOV, 
d = 100 yum, 



(11a) 
(lib) 



we find that go = 0.98. 

In addition to the trap strength, it is of practical inter- 
est to know the depth of the trap, that is, the minimal 
energy required for a trapped ion to escape. To establish 
the depth we assume that the amplitude of the micromo- 
tion induced by the rf field is small compared with the 
dimensions of the potential. In this c;asc, the depth of 
the trap region corresponds to the height of the (lowest) 
saddle point, Vg, of the trapping (local) minima of Ugs- 
We will for now ignore the control-field contributions to 
Ues, and only consider the "intrinsic depth," Us = Up{rs) 
provided by the ponderomotive potential of the rf field. 
By Eq. (1), Us is equal to 



Us= Uo 



Eriirs) 



Vri/d 



Uo = 



Q 



2y2 
rf 



4Mf22d2' 



(12) 



where Uq has the scaling properties of the depth. As 
we shall see by example in Sec. VB, control fields can 
significantly change the effective depth of SE traps. 

For the example parameters of Eq. (11), Uq = 6.1 eV, 
which is not a typical depth: as illustrated by Fig. 3 the 
maximal obtainable depth for a ring trap corresponds to 
less than 2% of Uq. 



A. Ring -trap characteristics 



IV. SE MULTIPOLES 



The different values of the parameter appearing in 
Eq. (8) lead to traps with different properties. In the 



Linear rf traps, such as the linear Paul trap, are based 
on an rf multipole field, which is translationally invariant 
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by 



FIG. 4: (Color online) The potential (a) and field (b) above 
a strip electrode at potential V (rectangle) with edges at ri 
and r2 in an infinite grounded plane, (a) Since the solid an- 
gle spanned by the strip is = 2A(j), it follows from Eq. (4) 
that the equipotential lines are circular arcs through the strip 
edges as discussed in the text, (b) By Eq. (5), the field is pro- 
portional to the magnetic field that would result from running 
a current along the edges of the strip in directions -E'(^) 
is given algebraically by Eq. (13). The field (potential) can 
be compactly expressed as a complex Polya field (complex- 
valued potential) according to Eqs. (16) and (17). 



along the trap axis thus providing a multipole guide. In 
addition to the ponderomotive potential provided by the 
rf field, quasistatic control fields are used to control the 
total trapping potential as illustrated in Fig. 1. With 3D 
electrode structures, the multipole guide can typically be 
designed by symmetry considerations alone. In the linear 
Paul trap, four parallel electrodes are used, an approach 
which can be extended to higher-order multipoles [14]. 
The same symmetry considerations are not applicable to 
surface-electrode traps where only a few configurations 
of high symmetry have been analyzed [4, 9]. 

In the following, we aim to present a simple geomet- 
ric theory of multipole guides in SE traps, as well as 
analytical methods for calculating the field in general 
translationally symmetric electrode configurations. The 
program will be similar to that followed in Sec. Ill for 
the ring trap: We will first introduce a constructive 
parametrization of all possible SE multipoles with a spe- 
cific center position, and subsequently analyze trapping 
strength and depth in terms of the design parameters. In 
this case we will also briefly discuss the effects of biasing 
the rf electrodes. 



E{r) 



r - r,; 



(13) 



for r in the trap region where q is z {—z) for the left 
(right) edge as illustrated in Fig. 4(b). In terms of the 
azimuthal angle Aip spanned by the strip electrode as 
seen from r, the sohd angle spanned by the electrode is 
n = 2A0 so that by Eq. (4) V{r) = VA^/tt. This, in 
particular, implies that the equipotential lines are circles 
through ri and r2 as illustrated in Fig. 4(a). 

Although compact, the vector nature of Eq. (13) leads 
to somewhat unwieldy expressions for the combined field 
of several strip electrodes. A simpler representation can 
be obtained by mapping the point r in the x-y plane to 
the complex point p, 

p{r) = x + iy, (14) 
and the vector field E{r) to the complex function 



(15) 



As discussed in more detail in Appendix B, £J is the Polya 
field of E* and E* (p) is an analytical function of p if and 
only if £J is a divergence-free and irrotational vector field. 

In terms of E, the field of a strip electrode as given by 
Eq. (13) is 



E*{z) 



ITT 



Z- Zi 



(16) 



where z — p{r), Zi — p{ri), and qi = ■ z. We have 
expressed E* in terms of z rather than p because it turns 
out that Eq. (16) as well as Eq. (17) below hold with- 
out any modification for the boundary conditions of a 
grounded cylinder (z — c) as we shall see in the next 
section. 

The complex potential <i>(z) fulfilling E*{z) = 
—d^{z)/dz is given by 



$(z) = V-^qiln{z- z,), 

TT ^ ^ 



(17) 



A. Polya field 

In the remainder of the paper, we will orient the coor- 
dinate system so that the translational symmetry is along 
the z axis and the electrode plane is parallel to the y axis 
at X = d. In figures we will orient the coordinate system 
so that the origin (trap center) is above the electrode 
plane, i.e., with the x axis pointing downward. 

It follows from Eq. (5) that the field of a strip electrode 
with edges intersecting the x-y plane in ri and r2 is given 



up to a complex constant. The real- valued physical po- 
tential V is equal to the real part of <i> as discussed in 
Appendix B. For the purpose of calculating derivatives 
Eq. (17) is valid as it stands, but care should be taken to 
choose the branch cut for the logarithm to ensure conti- 
nuity if the value of <i> is to be calculated directly. 



B. Parametrization of SE multipoles 

With the Polya field description of strip electrode fields 
provided by Eqs. (16) and (17), identifying the field ze- 
ros, i.e., trapping points, for a set of n strip electrodes is 
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FIG. 5: (Color online) A graphical representation of the con- 
formal map c — > p{c) mapping the inside of a conducting 
cylinder to the trapping region above an infinite conducting 
plane, as defined formally by Eq. (18). Green circles and 
triangles indicate the fixpoints at the point on the electrode 
plane right below the trap center {c = p = d) and at the trap 
center (c = p = 0), respectively. Thick blue lines indicate 
the electrode surface in the cylinder picture (a) and the SE 
picture (b). 




Re 



FIG. 6: (Color online) Parametrization of multipole elec- 
trodes and geometrical implementation of the map Eq. (18) 
for the electrode surface. The plot shows the electrodes in the 
c and p representation superimposed. Electrode positions for 
SE multipoles are specified by the parameters n, 6o, and 0^ 
according to Eq. (20). The example configuration corresponds 
to that of Fig. 7(b). An arbitrary point a = dexp{i(fii) on the 
electrode plane of the cylinder is mapped to pi = p{ci) on the 
SE electrode plane Re(p) = d according to Eqs. (18) or (19). 
The mapping can be performed geometrically by extending 
the circle tangent at Ci to the line Re(p) — d (dash-dotted), 
or by a Mercator projection to the imaginary axis (dotted). 



ations alone as detailed below. 

We choose the trap center as the origin at a distance 
d from the electrode plane given hy x = d, or equiva- 
lently, Re(p) = d, as illustrated in Fig. 5. A conformal 
map which maps the cylinder |c| — d to the electrode 
plane Re(p) = d with fixpoints at c=p=0 (trap center) 
and c=p=d (point on the electrode plane below the trap 
center) is given by the conformal (Mobius) map 



as illustrated in Fig. 5. 

For the purpose of trap design, we are mostly con- 
cerned with the mapping of the electrode structure. Here, 
we find that a point c = cZexp(i(/)) on the cylinder surface 
is mapped to 

r(p(de^'^)) =da; + dtan(f)y, (19) 

with a simple geometrical interpretation as illustrated in 
Fig. 6. 

For the purpose of characterizing the SE multipoles, 
we will introduce a parametrization of the simplest pos- 
sible multipole guides where an nth-order guide is im- 
plemented with n strip electrodes placed equidistantly 
around the cylinder walls. We will parametrize this con- 
figuration by the angle of 9^ spanned by each strip as 
seen from the cylinder center, and the azimuthal angle 
00 of the center of one strip, so that the electrode edge 
positions are at c = dexp{i(j>i^'^) for ^l^'' given by 

</>i+) = 0o- V + f {l,2,...,n}, (20a) 
(/.^) = 0o+ V + f {l,2,...,n}, (20b) 

as illustrated by Fig. 6. 

Together, Eqs. (19) and (20) give the explicit position 
of the edges of strip electrodes to create an nth-order 
SE multipole guide at a distance d above the surface as 
illustrated in Fig. 7. 

C. Field of SE multipoles 



reduced to the problem of finding the roots of a polyno- 
mial of order 2n — 1 . Similarly, the task of characterizing 
the identified trapping points in terms of multipole co- 
efficients is reduced to series expansion as discussed in 
Appendix B. 

Although this can be achieved in many cases, we are 
often interested in the inverse problem: how to place strip 
electrodes to achieve a multipole guide of certain strength 
and orientation, and how these parameters relate to total 
electrode area and other implementation considerations. 
To solve this, it is convenient to map the trapping re- 
gion conformally to a cylinder centered at the desired 
guide center, as an arbitrary multipole in a cylinder with 
conducting walls can be designed by symmetry consider- 



To determine the analytical form of the field in an SE 
multipole guide described by Eq. (20), we will calculate 
the field in the corresponding cylinder configuration and 
then calculate the SE field through the map Eq. (18). 

Firstly, we will determine the field of a strip electrode 
with edges at ci, C2 in a cylinder with conducting walls. 
Calculating i>(p(c)) with pi — p{ci) by Eq. (17) through 
the map Eq. (18), we find that when considering pairs 
of edges with opposite sign qi + q2 — 0, <i>(c) is equiv- 
alent to the value obtained directly from Eq. (17) with 
z — c and Zi — Ci, as illustrated by Fig. 8. In other 
words, Eqs. (13), (16), and (17) for the field and po- 
tential of strip electrodes in an infinite grounded plane 
hold without any modification for the case of similar strip 
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FIG. 7: (Color online) Examples of SE multipoles. Rectangles indicate rf electrode positions. Solid lines are contour lines of 
Up, spaced by a factor of \/2. The outermost line is at twice the maximal obtainable depth for quadrupole traps Eq. (30), 
corresponding to 2 x 55.8 meV for the example parameters of Eq. (11). Dashed lines are field lines (some omitted). Insets 
show the electrode configuration in the cylinder representation. Triangle and star marks indicate the guide center and intrinsic 
saddle, respectively. In terms of the parametrization Eq. (20), the configurations are described by (a) the symmetric quadrupole 
(n = 2) configuration (^o = 'r/2) with 6^ chosen for optimal intrinsic depth according to Eq. (28); (b) a tilted quadrupole 
configuration with 9^ — 7r/4 and = 57r/8; and (c) an octupole (n = 3) configuration with 6^ — tt/8 and 6o = n/16. 




FIG. 8: (Color online) The electric field of the field of a strip 
electrode at potential y in a conducting cylinder has the exact 
same algebraic form as the field of a strip electrode in an 
infinite grounded plane, as given by Eq. (5). 



electrodes on the walls of a grounded cylinder. This is in 
agreement with the fact that the field of an SE strip elec- 
trode meets the boundary conditions for a strip electrode 
in a conducting cylinder as illustrated by Fig. 4. 

We can write the values of the electrode 

edges Eq. (20) as 0^°^ + 9o T where = 

2tt{1,2, . . . ,n}/n. For {ci} — dexjp{i(j}n^) and q-i — 1, 
we find that the sum Eq. (17) can be completed since 
11^=1 i'^ ~ e'^'^'^/") = a" — 1, and the corresponding (un- 
physical) single-sided complex potential is given by 



<°nc) = V^^ln((i)"-l) 



can 



Since the positive and negative edge positions 
be obtained from ^i^-* by rotations through an angle of 
9q T (^w around the origin, the combined potential <i>„(c) 



is equal to ^l*^' (ce 



./2) 



ce 



-iieo+e^/2)) 



which we write as 



-i(eo-8™/2). 



(21) 



Equation (21) gives the exact form of a the complex 
potential of a cylinder multipole as parametrized by 
Eq. (20), and so together with the conformal map 
Eq. (18) give the potential for an SE multipole in the gap- 
less plane approximation. The physical potential with 
the specified boundary conditions is Re{<i>„[c(p(r))]} + 

To learn about the strength and orientation of the cor- 
responding multipole field, we expand 9<i>„(c)/9c to low- 
est order about the origin. 



dc 



V 2n 



■ sm 



c"-i +0(02"-!). (22) 



Since by Eq. (18) dc/dp = 1/2 at the origin, we have by 
Eq. (B7) that the lowest nonzero multipole coefficient of 
is 



2"" - sin 

TT 



n6^ 



V 



(23) 



with ap"'' as defined in Eq. (B6). The strength \a'^^\ is 
seen to reach its maximal value of 2V/7r(2c?)" for 9^ — 
7r/n, corresponding to half the cylinder being covered 
with electrodes. The complex phase factor exp(— m^g) 
describes the orientation of the multipole. 

At this point we can compare the strength of an SE 
multipole with that of a traditional (3D) multipole guide 
with the same ion-electrode distance. The strongest such 
guide is obtained when the electrodes are chosen as the 
ztV/2 equipotential surfaces of Re(z"V/2(i"), correspond- 
ing to a multipole coefficient of 



(") 

^3D 



1 V 

2~(f' 



(24) 
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so that the strongest possible SE multipole is a factor of 
2"7r/4 weaker than the strongest possible 3D multipole of 
the same order and with the same ion-electrode distance. 

For the case of quadrupole (n — 2) guides em- 
ployed in Paul traps, the strength of the rf field is 
usually parametrized by the dimensionless parameter q 
[10, 12, 14], which is related to the multipole expansion 
coefficients by 



9 ^ 



v(2) 



4 Q 



2tt 



10, 



(25) 



for qo defined by Eq. (9). In the limit where the adia- 
batic approximation is valid, the corresponding secular 
frequency (in the absence of bias fields as discussed in 
Sec. VB), is well approximated by q^/V& [12], corre- 
sponding to a maximal secular frequency of 5.5 MHz for 
the example parameters of Eq. (11). To compare the re- 
sults presented here with experimental results, we note 
that Ref. [5] gives a measured value of ^exp = 0.54 in the 
adiabatic approximation for an SE trap with 9^, = 7r/4, 
6*0 = 7r/4, M = 24 AMU, V,{ = 103 V, and d = 40 /im. 



This is in good agreement with the value q 
dieted by Eq. (25). 



0.55 pre- 




FIG. 9; (Color online) Intrinsic depth of quadrupoles (n = 2) 
as a function of 6o for Ow/ir equal to 0.3, 0.4, and 0.5 with 
the parametrization described by Eq. (20) for the example pa- 
rameters of Eq. (11). The depth scales as Uo (Eq. (12)). Solid 
lines are exact results, corresponding to cuts of the contour 
plot shown in Fig. 10. Dashed lines correspond to evaluating 
the exact field at the simplest saddle-point estimate pi'^'"'' of 
Eq. (26). Higher-order approximation as given by Eq. (A3a), 
could not be discerned from the exact results. Dash-dotted 
lines show the value of the crude estimate of Eq. (27) . 



V. DEPTH OF SE MULTIPOLE TRAPS 

In this section we will investigate the potential depth 
that can be obtained when using SE multipoles for rf 
traps. As in Sec. Ill on the SE ring trap, we will as- 
sume the adiabatic approximation Eq. (2) to be valid. 
In Sec. V A, we will consider the intrinsic depth, i.e., the 
depth in the absence of control fields. Here we will mostly 
work with crude estimates, while Appendix A provides a 
derivation of some exact results. In Sec. VB we consider 
the effects of control fields on the total potential depth 
provided by Ucs- 

Our current understanding of SE multipole depth is 
far from complete and, in particular. Sec. VA contains 
results based on conjectures as detailed in Appendix A. 



A. Intrinsic depth 

In terms of the complex rf potential <i>i.f the pondero- 
motive potential is proportional to | [^ where we use 
a prime to denote differentiation with respect to p. It 
follows that a necessary condition for an intrinsic saddle 
point is that <&J.f'&"f = 0- Since the zeros of <l>J.f are not 
saddles, the saddle point ps must fulfill <i>"f(ps) = 0. 

For the multipole configuration described by the com- 
plex potential <!>„, we can take advantage of the fact that 
the expansion <!>„ « ac"''c" by Eq. (21) is good to order 
2n — 1 to establish a reasonable estimate of the saddle 



point pi"-* fulffiling <i>"(ps"-') — 0. By the chain rule 



(")^ _ 



^„ _ 6)2$ _ a$ d^p I dp 

dp'^ \ dc^ dc dc^ I dc 



so that for <i> cx c" we have a saddle at c = d{l—n) /{1+n) 
corresponding to p = d{l — n), which we take to be a 

reasonable estimate of pi"'' as follows: 



P 



,(«) 



,(»,o) ^ 



= d{-n+l). 



(26) 



The saddle point of Up for is consequently located ap- 
proximately above the center of the guide at a distance 
of n times the ion-electrode distance from the electrode 
plane. An iterative procedure for obtaining better esti- 



mates p'f''^\ k 



1,2,... of the true saddle Ps' is given 
by Eq. (A3a). 

As illustrated by Fig. 9, Up ^pi"'"''^ is a reasonable esti- 
mate of the intrinsic depth [/i"^ = Up {ps^^^ correspond- 
ing to While Up can be evaluated exactly according 
to Eq. (21) by Eq. (Al), the expression is somewhat in- 
volved. A crude estimate of Us^'^ can be obtained by 
employing the approximation <&,i(p) ~ ac"^c(p)" again 
to estimate the <i>'(pi"'"-'), resulting in 



V/d 



(27) 

which for the example parameters of Eq. (11) corresponds 
■sin^ (n 6*^/2) 180 meV. The value of the 



to U, 

crude estimate is illustrated in Fig. 9. 

A more detailed analysis of the intrinsic depth is car- 
ried out in Appendix A. The main conclusion is that for 
any value of n there is a special saddle position pi"'' and 
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FIG. 10: (Color online) Contour plot showing the intrinsic 
depth Ub"^ of quadrupole (n — 2) guides as a function of the 
parameters and 6^ appearing in Eq. (20). Contour lines are 
spaced by 0.1 times the maximal quadrupole depth Us given 
by Eq. (30). The dashed line indicates points where Eq. (28) 
holds. For these points the maximal depth Us"^ = [/i"' is 
obtained. The symmetries of the depth are discussed in Ap- 
pendix A, and the dash-dotted lines (dots) indicate parameter 
values also illustrated in Fig. 9 [Figs. 7(a) and 7(b)]. 



a real-valued constant A. 
rameters 9o and 9- 
and Oo fulfill 



both independent of the pa- 
if and only if 9^, 



so that pi"-* ~ pi"'' 



cos{n9w/2) ~ AnCos{n9()). 



(28) 



Further, the intrinsic depth is the same for all com- 
binations of 6*0 and 9yj fulfilling this requirement, and 
this depth is the maximum obtainable in an nth-order 

~ (n) 

SE multipole. To obtain the maximal depth Us = 
(pi"^) and maximal strength simultaneously, we 
must, by Eq. (23), have 9^ = T^/n, so that according 
to Eq. (28), cos(n^o) must also vanish. This corresponds 
to an antisymmetrical configuration where a refiection in 
the real axis (z axis) map rf electrodes to ground elec- 
trodes and vice versa. 

For quadrupoles (n = 2) , A2 — VE — 2 and the special 
saddle is located at 



P. 



(?) ^ -d 



/5-I 



(29) 



The maximal intrinsic depth of a quadrupole trap is 

11 



27r2 



-f/o, 



(30) 



with Uq given by Eq. (12). For the example parameters 

of Eq. (11) Us — 55.8 meV. In the quadrupole case, 
the antisymmetric configuration (9q = 7r/2, 9^ = 7''/2) 
achieving maximal depth and strength simultaneously is 
referred to as the four- wire trap [4, 5]. 



B. Depth in the presence of static fields 

The remainder of the section is devoted to a brief de- 
scription of the impact on static control fields on trap 
depth. For simplicity, we will not consider the effect of a 
confining potential along the axis of translational symme- 
try. Also, we only consider fields that vanish at the trap 
center, since any field at this point would shift the trap- 
ping minimum of UcS to a region of nonzero rf amplitude 
leading to micromotion. 

By applying different control voltages to different elec- 
trodes, a control field of this type can be implemented 
with almost any set of two or more electrodes. Here, we 
will only consider the possibility of applying a bias volt- 
age V^crf to the rf electrodes. In practical applications, 
this approach has the advantage that a control field gen- 
erated by the rf electrodes is guaranteed to have a field 
zero exactly at the trap center, even when departures 
from the idealized gapless plane geometry are taken into 
account. Also, this approach ensures that the principal 
axes of the control and rf quadrupoles are aligned, al- 
lowing an analytical solution for the motion of a trapped 
ion valid outside the adiabatic approximation [14]. The 
main practical obstacle to applying a bias to the rf field 
is that it is incompatible with grounding part of the rf 
resonator. 

For a bias voltage Vai applied to the rf electrodes, the 
combined effective potential is equal to 



(31) 



where 6rf(T') = K-f('")/K-f is the rf electrode basis func- 
tion and 



c.rf 



(32) 



parametrizes the strength of the bias field. The ratio of 
the modified to the intrinsic depth only depends on Vc, 
so that the relative improvement of trap depth, which 
can be achieved by biasing, is a purely geometrical factor 
and does not depend on the operating parameters of the 
trap. 

An important consideration when applying a control 
field is the effect on trap stability. For the quadrupole 
trap with a control field applied as a bias to the rf 
electrode, the stability criteria can be established ana- 
lytically [10, 14]. For this purpose, the bias strength 
is parametrized by the dimensionless parameter a = 
2 9K;,rf/K-f together with q as defined by Eq. (25). The 
stable region usually used for QIP ion traps includes 
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a) ^ ^ b) 




FIG. 11: (Color online) Contour plot of UeH for the unbiased 
(a) and optimally biased (b) quadrupole configuration with 
9o — 100° and 6^ — -k 12. In order to illustrate the potential in 
the full trapping region, the potential has been mapped to the 
c plane so that what is plotted is f/cff (p(c)) with p(c) given by 
Eq. (18). Contours are spaced by a factor of \f2 from [/i'^Vl 
to 64[7i^-' with tjf' as defined in Eq. (30). Stars indicate the 
lowest saddle points, and triangles indicate the trap centers. 
As discussed in the text, the relative increase in depth by a 
factor of 9.8 is independent of operating parameters. 

(a, g) for which q < 0.7 and la/q^j < 0.5 [10, 14], where 
the second constraint is the one of interest here. If 
we introduce the dimensionless "geometrical" quadrupole 
strength a^^^ = a^^^'d'^/V, so that q — |a^^^|go, we have 
that a/g2 = wj8|a(2)| so that the stability requirement 
is a geometrical property independent of the operating 
parameters. For the parametrization of Eq. (20) we have 
by Eq. (23) that |a(2)| = sin(6i^)/27r so that the stability 
criterion is well approximated by |uc| ^ 0.6 sin (0^,). 

As an example, we consider a quadrupole guide (n = 2) 
with parameters 6*0 = 100°, and 0^ = tt/2. Here we 

- (2) 

find numerically that a maximum depth of 9. 8 Us is ob- 
tained for Vc — 0.18, corresponding to a/q^ = 0.14, well 
inside the stable region described above. The form of the 
optimally biased Ues is illustrated by Fig. 11. For the 
example parameters of Eq. (11), the optimal bias volt- 
age is Vc,rf = VcUq/Q = 1.1 V. Note that although the 
optimal bias voltage is dependent on operating parame- 
ters, the ratio a/q^ at optimal bias is constant, and so 
the optimally biased configuration will always be stable, 
provided it meets the constraint q < 0.7 corresponding 
to a lower bound on 17. 



VI. CONCLUSION AND OUTLOOK 

In conclusion, we have obtained a wide range of an- 
alytical results on SE traps. The most significant are 
the explicit parametrization of SE multipole geometry as 
given by Eq. (19), together with the results on the exact 
field [Eq. (21)], strength [Eq. (23)], and intrinsic depth 
[Eq. (30), Figs. 9 and 10]. 

There are, however, a number of open questions relat- 
ing to SE traps. 

In direct continuation of the work presented here it 



would be of interest to establish a formal proof that the 
intrinsic multipole depth is maximal exactly when the 
saddle is located on the special saddle point as conjec- 
tured in Appendix A. Also, the adiabatic approximation 
will in many cases be invalid within the stable region 
suggested by the intrinsic depth. It would consequently 
be of interest to study the depth beyond the adiabatic 
approximation [25]. 

The limitations of the gapless plane approximation 
used throughout this paper should be quantified in two 
respects. Least importantly, the effect of finite gaps be- 
tween electrodes are of interest because increasing the 
gap size decreases field gradients and the capacitive load 
of the trap. More important are finite-size effects: In a 
realistic implementation there will be grounded surfaces 
within a few cm of the electrode surface. Finite-size cor- 
rections will be first order in the total electrode extend 
relative to this distance, and could consequently be sig- 
nificant. 

Finally, limitations on the possible global structure of 
SE fields are not fully understood. In spite of recent 
progress on the possible structure of intersections for rf 
trap networks [11], it is not clear whether it is possible to 
implement ideal intersections in SE traps. Also, arrays 
of SE traps are being studied in several contexts [21, 26]. 
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APPENDIX A: SADDLES IN MULTIPOLE 
GUIDES 

This appendix is devoted to a detailed discussion of the 
exact intrinsic depth and saddle position in SE multipole 
guides. We use the term "intrinsic saddle point" to refer 
to a saddle point of the ponderomotive part of UcS as 
given by Eq. (2). As discussed in Sec. VA, such saddle 
points must be located at zeros of <i>"f. The following 
is devoted to identifying such zeros for <i>i-f given by the 
exact multipole potential <&„ of Eq. (1). To make the 
algebra more manageable, we will set V and d to unity, 
and introduce a new complex coordinate system with its 
origin in the electrode plane. 



With these conventions, — (d^n/dc) /{du/dc) can 
be given a compact form in terms of c and c — e~^^°c. 

^1 sin(n6>^/2) (1 + c)^ 

" TTc" -2cos(n6'^/2) + c-" c ' ^ ' 
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FIG. 12: (Color online) Orbits of the intrinsic saddle position 
for n = 2 and 3 as a So is varied for nQ^j-K = 0.7, 0.9, and 
1. Solid (dashed) lines show the exact value (the estimate 
ui"'^') of Us"'. Dots show the crude estimates ui"'"' = —n. 
Note that when is too large to be compatible with Eq. (A4) 
for any value of 6'o, the orbits do not pass through the special 
saddle positions ui"'' (stars). 



In terms of w, we have that c = (1 4- u)/(l — u) so that 



4n . / nBu 



P{u) 



where the polynomial P(u) = P+(w)-|-«P_(u) is described 
as a sum of even and odd parts, Pj-, given by 

P+(u) =cos(n0o) [(1 - + (1 + 

-2(l-u2)"cos(^^^ , (A2a) 

P_(?^) =sin(n6lo) [(1 - u)^" - (1 + u)^"] . (A2b) 

We see that <i>^ has zeros of order n— 1 at 7i = ±1, corre- 
sponding to the multipole center and its mirror image in 
the electrode plane. Also, there are 2n first-order poles 
at the roots of P, which are located at the electrode edges 
on the imaginary axis. 

To establish the saddle positions, we differentiate 
again, and write the derivative as 



(J)" — — <i)' 



1 



P(u)(l -1*2) 



where S(u) = S^(u) -|- iS-(u) is a polynomial of order 
2n -|- 1 with even or odd parts given by 

S±{u) = 2u(n - l)P±(u) -I- (1 - m2)P4(m), 

so that, in particular, 

S_{u) = -2sin(n6lo)(u {(u - 1)^" - (m + 1)^"] 

+ n[(?/~l)2" + (u+l)2«]). 

We see that <i>^ has zeros of order n — 2atu = ±l and 
second-order poles at the electrode edges, while the 2n-f 1 
roots of S describe the possible saddle positions. 



To determine the root structure for 5", we note the 
following: According to the edge-current picture, the field 
must be antisymmetric under inversion in the electrode 
plane so that <i>'(u)* = <&'(-«*). This implies that <!>' 
is real on the imaginary axis, or equivalently, that the 
electrode field is perpendicular to the electrodes at the 
surface as it should be. Since will diverge with equal 
sign at the two sides of an electrode or ground strip, 
there must be a zero of <i>" on each such segment on the 
imaginary axis. Further, since these zeros must also be 
roots of S{u) , there will be a zero of S between each of the 
2n edges for a total of 2n — 1 zeros on the electrode plane, 
i.e., the imaginary axis. As S(u) is a polynomial of order 
2n+ 1 this leaves two roots outside the electrode plane, 
which we identify as the saddle interest to us, which we 
will denote ui"'' and for which Re(Mi"'') < 0, and its 
mirror image —[Us ) ■ 

It follows from the above discussion that the saddle 
position ui"^ corresponding to <!>„ as given by Eq. (21) is 
the unique root of S{u) = with Re(u) < 0. The roots of 
S can be calculated numerically by standard techniques, 
allowing a simple procedure for determining the intrin- 
sic saddle points for any multipole. Alternatively, the 
following heuristic iterative procedure 



,(».o) _ 



1 ~ P'{u) 



(A3a) 
(A3b) 



found by solving S(u) = for the first order u, has been 
observed to converge relatively quickly to the main sad- 
dle as illustrated in Fig. 12. In particular, for estimating 
the trap depth, the zeroth- or first-order estimates es- 
tablished by evaluating <i>^ at ui^''^^ or ui^''^^ are mostly 
sufficient, as illustrated in Fig. 9. The initial estimate 
mI"'"'' is that identified in Sec. VA, but we note that it 
could also be obtained by starting the iteration at the 
trap center, u = —I. 



1. Saddle on the real axis 

It turns out that configurations where the saddle is on 
the real axis, i.e., directly above the trap center, have a 
special significance. 

There are two types of symmetry that will guarantee 
the saddle to be on the real axis. Firstly, when Oq is 
equal to an integer multiple of 7r/n, the electrode con- 
figuration is symmetric under refiection in the real axis, 
ensuring that the saddle must be located on the axis. In 
this symmetric configuration, S-{u) is identically zero so 
that S = S'_|_. A less obvious symmetry is found when 
dw — T^/n and 9q is an odd multiple of 7r/2n, so that the 
electrode structure is antisymmetric under refiection in 
the real axis, i.e., rf electrodes are mapped to grounded 
areas and vice versa. Since the ponderomotive poten- 
tial only depends on the field strength, this symmetry 
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TABLE I: Numerical values of the special saddle point wi"' 
and the factor appearing in Eq. (A4), describing the rela- 
tion b(>tw(X'n Oi) and 0,,,. for ojitiinal dejith. Exact values for 
n = 2 are given by Eq. (29). 



n 


_(n) , 
—Ua /n 


A 


2 


1.02909 


0.236068 


3 


1.03742 


-0.266149 


4 


1.04044 


0.276076 


10 


1.04375 


0.286475 


20 


1.04422 


0.287935 


50 


1.04436 


0.288342 


100 


1.04438 


0.288400 


200 


1.04438 


0.288415 



will also conserve the saddle which must consequently be 
located on the real axis. 

The antisymmetric configuration turns out to be of 
special interest, and we will refer to the saddle position 
in this case as wi"' . wi"' can be calculated as the negative 
real root of S- , since in the antisymmetric configuration 
cos{n9o) and cos{n0yj/2) both vanish so that S = iS-. 

Although we have introduced wl"^ as the saddle posi- 
tion for the anti-symmetric configuration, so that in par- 
ticular ul"-* is independent of and 9yj, it turns out that 
there is a number of configurations with the same saddle 
position, and that these have very interesting properties. 
Firstly, we note that since P+ and thus depends on 0q 
and Oy, only through the ratio of cos(n^,„/2) to cos(n^o)5 
there will for any value of n be a constant A„ so that 
5+ = whenever 



cos(n(?^/2) = ^„cos(n(?o) 



(A4) 



Since the roots of S- are independent of 9o and 9^, it fol- 
lows that when the condition of Eq. (A4) is met, the spe- 
cial saddle ui"-* will be a root both of 5+ and S- . Since 

S = iS^ it follows that ui"' will also be a root of S 
and consequently be the actual saddle position through 
an "accidental" symmetry as illustrated in Fig. 12. We 
have not been able to obtain closed expressions for ui"^ 
for n > 2, but give numerical values in Table I. 

The key property of the special saddle point is that 
all configurations which fulfill Eq. (A4), and thus have 



the saddle at 



obtain the maximal intrinsic depth 



possible for a SE multipole of order n. 

Ui^\e^,9o) = tJi''^ ^ ui"\9^,9o) = ui^l 

This result has only been proven for n = 2 but we con- 
jecture it to be true for all n based on numerical results. 

APPENDIX B: 2D ELECTROSTATICS 

This section recapitulates a number of results on 2D 
electrostatics. 



1. Polya field 

The Polya field gives a one-to-one mapping from 2D 

electrostatic fields to analytical functions on the complex 
plane [27] . We consider the canonical mapping (a;, y) i-^ 
z = X + iy from to C and map the electrical field 
E{x, y) to the complex-valued function 



E{x + iy) = E^{x,y) + iEy{x,y). 



(Bl) 



Note that E is not related to the complex amplitude used 
in electrodynamics to describe phase relations. 

We will consider z = x + iy and z = x — iy to he 
independent variables with respect to differentiation, so 
that = \dx — i\dy and dz = \dx + i\dy are the 
Wirtinger derivatives often used in physics. With this 
convention we find that 



dE 
dz 



dE* 
dz 



1 



{V ■ E + iz-V X E), 



(B2) 



so that V • E and 'V X E vanish, as they must for an 
electrostatic field in free space, if and only if dE/dz = 0. 
In terms of the Wirtinger derivatives, a complex func- 
tion /(z) is differentiable if and only if df{z)/dz = 0, 
and Eq. (B2) consequently implies that E is a. free space 
electrostatic field if and only if E* is complex differen- 
tiable and thus analytic. In complex analysis, the field 
E is known as the Polya field of E* . 

In terms of the real-valued potential V the analytical 
function E* is seen to be given by 



E* 



fdV 
\'d^ 



.dV 
dy 



^dV 
' dz ' 



(B3) 



Note that since V is real it is not analytic and espe- 
cially not an antiderivative of —E*/2. On the other hand, 
the analytic function —E* has an analytic antiderivative 
^{z),so that 



d^ d^ 

dz dz 



0. 



(B4) 



The antiderivative is unique up to the addition of a com- 
plex constant. By inspection, we see that V defined as 
V{x,y) = Rei^x + iy)) does fulfill Eq. (B3). The imagi- 
nary part of $ is known as the harmonic conjugate of V, 
and we note that since V Im(<i>) ■ E = 0, lines of constant 
Im($) are field lines of E. 



2. Conformal maps 

In two dimensions, conformal, or angle-preserving, 

maps will map a solution of the Laplace equation to an- 
other solution of the Laplace equation in the sense that 
if V'^V{r) = and M(r) is a conformal map, then we 
find that V^y(M(r)) = 0. In the complex plane a map 
w z{w) is conformal if and only if the function z{w) is 



13 



analytic. If E{z) is a divergence-free and irrotational 
field, so that E*{z) = —d^{z)/dz for some complex- 
valued analytic potential ^{z) and z = z{w) is an an- 
alytical function, ^^^^{w) = ^{z{w)) is also analytic so 
that by Eq. (B2) the Polya field of 



p*. X d<i>w{w) ~ ..dz{w) 
Ey,{w) = = E {z{w))- 



dw 



dw 



(B5) 



is indeed divergence-free and irrotational by the results 
of the last section. 



3. 2D multipole expansion 

In two dimensions, a naive specification of the nth. mul- 
tipole moment of a harmonic potential V would involve 
the n + 1 distinct partial derivatives of order n. Consid- 
ering the constraints imposed by the Laplace condition 
would, however, leave only two free parameters. Work- 
ing in the complex notation, this result is obvious, as the 



multipole expansion of V is nothing more than the series 
expansion of a corresponding analytical potential <1> so 
that V — Re(<i>). We will write the multipole expansion 
around the origin as 



a 



(B6) 



n=0 



We will now consider the mapping of the multipole 
strengths: If $(z) = + 0(2"+!) and = 

<i>(z(?x;)), where z — z{w) is a conformal map so that 
z(0) = 0, we see that <^^{w) = a^'^\dz/dwYw'' + 
0(11;"+^) so that the lowest-order multipole moment of 



dzY 

dw I 



(B7) 



where dz/dw is evaluated at the origin. 
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